Break-up of suspension drops settling under gravity in a viscous fluid close to a vertical wall 



O 

(N 



>v 
i 

q, 

o ■ 

• i-H , 
^ ■ 

Oh 1 



^1" 

o 
in 

(N 

o 



Anna Mytyk*, Walter Meile^, Giinter Brenn^ and Maria L. Ekiel-Jezewska* 
* Institute of Fundamental Technological Research, Polish 
Academy of Sciences, Pawihskiego 5B, 02-106 Warsaw, Poland 
^Institute of Fluid Mechanics and Heat Transfer, University of Technology, Inffeldgasse 25/F, 8010 Graz, Austria 

(Dated: January 13, 2013) 

The evolution of suspension drops sedimenting under gravity in a viscous fluid close to a vertical 
wall was studied experimentally and numerically with the use of the point-force model, in the Stokes 
flow regime. The fluid inside and outside the drop was identical. The initial distribution of the sus- 
pended solid heavy particles was uniform inside a spherical volume. In the experiments and in the 
simulations, the suspension drops evolved qualitatively in the same way as in an unbounded fluid. 
However, it was observed, both experimentally and numerically that, on the average, the destabiliza- 
tion time T and the distance L traveled by the drop until break-up were smaller for a closer distance 
h of the drop center from the wall, with approximately linear dependence of T and L on D/h, for h 
larger or comparable to the drop diameter D. Destabilization times and lengths of individual drops 
with different random configurations of the particles were shown to differ significantly from each 
other, owing to the chaotic nature of the particle dynamics. 



I. INTRODUCTION 

Motion, deformation and break-up of suspension 
drops settling under gravity in a viscous fluid has 
been observed and studied for many years by differ- 
ent researchers. 1-7 The most interesting feature of such 
a process is that a suspension drop settling under the 
influence of gravity in an unbounded viscous fluid re- 
mains a cohesive entity for a long time, even in the ab- 
sence of a surface tension at the drop surface, and with- 
out any attractive direct interactions between the sus- 
pended particles. 3,5 This property is interesting as a chal- 
lenging basic theoretical problem, as well as a typical 
process observed in many practical contexts, such as de- 
position of drugs in the human airways, 8 rising of mantle 
plumes, 9 or interaction of melamine formaldehyde with 
a complex plasma^ 

There is a lot of interest in studying deformation and 
destabilization of suspension micro-droplets, especially 
such which consist of a cloud of non-Brownian heavy 
particles which are denser than the fluid and are sepa- 
rated from each other by a fluid identical to the host fluid 
outside the drop^"— For such systems (on which we focus 
in this paper), the fluid inertia is negligible, and the fluid 
flow can be described by the Stokes equations. 11,12 Based 
on the similarity lawr^ the dynamics of such micro- 
systems is identical to the dynamics of geometrically 
similar, but larger and faster objects which move in a 
fluid of a higher viscosity, provided that the Reynolds 
number Re for both systems is the same. Suspension 
drops settling in very viscous fluids at Re < < 1 have been 
studied experimentally by standard video tracking^ 

The typical evolution pattern of a single suspension 
drop settling far from walls or interfaces is the follow- 
ing. Initially, the injected drop tends to become spher- 
ical, loosing a large number of particles, which move 
more slowly and therefore form a thin tail behind (above) 
the drop. The particles suspended in the drop recircu- 



late, settling down faster than the center of mass in the 
inner part of the drop, and then moving to the outer part 
of the drop and settling there more slowly than the cen- 
ter of mass, then coming back to the inner part, and so 
on. At the top of the drop, some of the particles sep- 
arate out; they move more slowly than the drop and 
form a thin tail above it. The drop slowly expands side- 
ways and becomes more and more flat. A hole inside is 
formed, and the resulting torus grows horizontally and 
decreases its height. Then, the drop suddenly breaks 
into two (or sometimes more) smaller droplets, which 
repeat the same evolution pattern. 

This scenario is qualitatively reproduced with the 
use of various numerical techniques for solving the 
low-Reynolds-number hydrodynamic equations: the 
simple point-particle model?' 5 or the Rotne-Prager 
approximation 6 - (for Re = 0), and the lattice-Boltzmann 
method, 1 * spectral methods^ 5 - or Oseen interactions^ (for 
a finite Reynolds number). 

The goal of the present paper is to investigate experi- 
mentally and theoretically the influence of a solid verti- 
cal wall on the evolution, deformation and destabiliza- 
tion of the suspension drops described above, with the 
emphasis on the destabilization time and the distance 
a suspension drop travels before breaking (destabiliza- 
tion length). The motivation comes from the previous 
work, 16 where simple experiments confirmed that the 
presence of one or two parallel vertical walls reduces the 
destabilization length. The effect of the wall on the desta- 
bilization of suspension drops is important for practical 
applications. 

The point-particle model in an unbounded fluid was 
successfully applied in the literature to determine basic 
features of a suspension drop evolution^"— Therefore, in 
this work, we use as the theoretical framework the point- 
particle model close to a plane solid wall, parallel to the 
direction of gravity. This model has been recently devel- 
oped and tested elsewhere. 17 To perform the measure- 
ments, we use the experimental set-up applied in a num- 



ber of previous publications^^ 

In Sec. HH we describe our experimental setup. In 
Sec. Hm it is used to generate data on the evolution pat- 
tern of suspension drops and their destabilization times 
and lengths. In Sec. IIV1 we introduce a simplified geo- 
metrical model of the internal structure of the suspen- 
sion drop, to allow for efficient numerical modeling. 
In Sec. |V] we introduce the point-particle model close 
to a wall, present numerical simulations of suspension 
drops, and determine destabilization times and lengths. 
Secs. lVII and lVIII are devoted to comparison between the 
measurements and the theory, discussion and conclu- 
sions. 



II. EXPERIMENTAL SET-UP 

Experiments were performed in a glass-wall ves- 
sel with the square cross-section 20 cm x 20 cm 
and the height 100 cm (see Fig. [TJ. The container 
was filled almost up to the top with 86% glycerol 
of density p = 1224 kg/m 3 and kinematic viscosity 
v = 0.9-10 -4 m 2 / s. The corresponding dynamic viscosity 
?/ = pv = 110 mPa s. 

A suspension was prepared by mixing glass particles 
with the 86% glycerol (the fluid taken from the con- 
tainer). The volume fraction was set constant to <f> = 0.1. 
To control the volume fraction, a certain quantity of glyc- 
erol was weighted to determine its mass M. The mass 
Mp of glass particles required to obtain the volume frac- 
tion <f> = 0.1 was calculated from the relation 

_ Mp/pp 
9 M/ p+ Mp/pp 1 K ' 

with the glass-particle density pp = 2400 kg/m 3 . 

The glass particles were polydisperse. The distribu- 
tion of their radii was measured using a Nicon Eclipse 
E50 optical microscope equipped with a 20 x magnify- 
ing lens. A typical image observed under the micro- 
scope is shown in Fig. [2] The actual size of 411 particles 
was measured, resulting in the mean particle diameter 
d = (25 ± 8)/jm. 

A drop of the suspension was generated by a trigger 
mechanism, shown in the top panel of Fig.Q] This mech- 
anism consisted of a pump-needle-unit (enlarged in the 
bottom panel of Fig. [I}, a sequencer, a function generator 
and a power supply. The sequencer allowed to control 
the size of droplets. In the experiments, the drop diam- 
eter was close to 1 mm. 

A single suspension drop was released at a distance h 
from the side wall of the container and at the distance 
z a = 4 mm below the free surface of the fluid. Its mo- 
tion and shape evolution were observed, as illustrated in 
Fig. |3 and recorded by a SONY XCD-X710 video cam- 
era (black and white). The position of the camera was 
changed stepwise to keep the drop within the field of 
view. 




FIG. 1: The experimental set-up: the drop generator with its 
driver, the container, and the devices for visualization (top). 
Photograph of the drop generator (bottom). 

The goal of our investigations was to study the influ- 
ence of a vertical wall on the motion, shape and break- 
up of a suspension drop. Therefore, we investigated the 
evolution of drops, which were released at different dis- 
tances h = 2,5, 10, 20, 30, 50, 100 mm between their cen- 
ters and the wall. For a given value of h, the measure- 
ments were typically repeated n = 20 — 35 times, to de- 
termine the dispersion of the results and reduce statisti- 
cal errors. In total, 202 measurements were performed. 

To control the initial size of all the drops, at each exper- 
iment the procedure of release was tested by measuring 
the time t the drop needed to settle the same distance 
y, (with y = 30 cm for h = 2 mm and h = 5 mm, and 
y = 50 mm for all other distances h from the wall), start- 
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FIG. 2: Glass particles under the optical microscope. 
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FIG. 3: Illustration of the experiments. 



ing at {z\ + z ) = 41 mm from the free surface of the fluid 
in the container, see Fig. |3] The actual initial diameter 
D of each drop was later determined from photographs 
taken with the same SONY video camera during the ex- 
periments, with the average over all the measurements 
(D) = 1.13 mm. 

The Reynolds number based on (D) and the velocity 
v = 0.92 mm/ s of the drop in the central part of the con- 
tainer was much smaller than unity, Re = (D)v/v « 10~ 2 . 

The number of particles iVo inside a single suspension 
drop of diameter D was estimated in the following way. 
First, the sizes of rid = 411 glass particles were measured 
under the microscope, and the total volume Vd of these 
rid particles was calculated. Since the volume Vd of all 
the No particles inside a suspension drop of diameter D 



and the volume fraction <f> = 0.1 was known, 

nD 3 



V D 



(> 



(2) 



then for each suspension drop, No was determined from 
the rescaling relation, 



N 



Vpn d 
V d 



(3) 



Finally, the number of particles in each suspension 
drop Nq was averaged over all the experiments, with 
(N ) = 7100 ± 200. 



III. EXPERIMENTAL RESULTS 
A. General observations 

During the experiments it was observed that a suspen- 
sion drop settling close to a vertical wall evolved qualita- 
tively in a similar way as in an unbounded fluid. At the 
beginning, the drop remained almost spherical, leaving 
a thin tail of particles behind. The drop shape gradu- 
ally changed into a horizontal torus, which continued to 
flatten and widen, and suddenly deformed and destabi- 
lized into two (or sometimes more) droplets. Therefore, 
the main characteristic features of the drop evolution ob- 
served previously in an unbounded fluid (see, e.g., Fig. 
23 in Ref. H and Figs. 18 and 21 in Ref. H) are also 
seen in the presence of a vertical wall. There are, how- 
ever, the following differences. Close to a wall oriented 
parallel to the gravitational field, the drop looses more 
particles; also, it flattens and breaks up earlier. More- 
over, the drop settling very close to the wall (at h = 2 
mm) destabilized in such a way that the line of centers 
of the two droplets formed after the break-up in most 
cases was oriented along the wall (in 23 out of 25 such 
experiments). At ft. > 5 mm, the lines of centers were 
oriented randomly. 

The statistics of the observed evolution patterns is 
listed in Table|U All the drops broke up, 99% of them into 
two fragments, and only two drops into three pieces. 



Dependence of destabilization time and length on the 
distance from the wall 



In the experiments, we measured the instant T of drop 
destabilization, and the vertical drop position L at that 
time, as functions of the distance h from the wall. The 
destabilization time T was defined as the moment when 
the torus began to break, and the destabilization length L 
as the distance traveled by the drop from its release (cor- 
responding to time equal to zero) until destabilization, 
see Fig. [3] Breaking of the torus was observed on the 
side view. Just before destabilization, the torus became 
thinner in the central part than at the outer ones. Then, 
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TABLE I: The number of suspension drops with the indicated 
final stage of the evolution in our experiments. 



h/D 


destabilization 


destabilization 


no break-up 


total 




into 2 droplets 


into 3 droplets 






87.7 


34 








34 


43.1 


35 


1 





36 


26.8 


33 








33 


17.9 


27 


1 





28 


8.6 


22 








22 


4.6 


21 








21 


1.8 


24 








24 



it began to bend, with the thinner central part above 
the thicker outer ones, and this moment was called the 
destabilization time. For each drop, T and L were mea- 
sured using a stopwatch and a scale on the back wall of 
the container. Then, values of T and L were averaged 
over all n experiments performed for a given distance h. 

Values of the average destabilization time and length 
for different distances h are listed in TablelTTl For simplic- 
ity, we use here the same symbols T and L, which from 
now on in this paper will denote the average values. 



TABLE II: The destabilization time T and length L of suspen- 
sion drops as functions of the distance from the wall h, aver- 
aged over all experiments performed for a given distance h. 
The average initial drop diameter D and the inverse time unit, 
1 /t, defined in Eq. (6j, are also indicated. 



h [mm] 


D [mm] 


L [mm] 


T[s] 


1/r [1/s] 


2 


1.14 


186 ± 10 


348 ± 19 


1.59 


5 


1.08 


266 ±9 


434 ± 16 


1.51 


10 


1.16 


326 ±8 


501 ± 10 


1.62 


20 


1.12 


319 ± 10 


477 ± 13 


1.56 


30 


1.12 


333 ±7 


467 ± 11 


1.57 


50 


1.16 


344 ±8 


466 ±9 


1.61 


100 


1.14 


342 ±8 


488 ± 10 


1.59 



drop initial diameter D, its initial settling velocity 

v — ^ tot 



(4) 



determined by the total gravitational force (weight mi- 
nus buoyancy) acting on the drop, 



Ft. 



-(pp - p)g, 



and the corresponding time unit, 

D 

T = 2V 



(5) 
(6) 

Here g is the 



or explicitly, r = 15r][2g(pp — p)D<p]^ 1 
gravitational acceleration. 

In Figs. H]and|5j the average non-dimensional destabi- 
lization length L/D and time T/t are plotted versus the 
non-dimensional distance h/D of the drop center from 
the wall. In both figures, the statistical errors are indi- 
cated by the black error bars. The measurements of time 
T/t have an additional "systematic" error, related to un- 
certainty of t, owing to fluctuations of the initial drop ve- 
locity v from day to day. They might be related to small 
temperature variations of the fluid, temperature gradi- 
ents, convection of the fluid, water at the interface of the 
glycerol, or other reasons. The method to estimate such 
a systematic error, shown in Fig.|5]by the gray error bars, 
is described in Appendix lAl 
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The actual initial diameter D of each drop was deter- 
mined from photographs taken during the experiment 
at a time instant between 3 s and 12 s, and averaged over 
the experiments performed at a given distance h. If the 
diameter of a drop was more than 15 percent different 
from the average, such an experiment was discarded. 
Then, the average diameter was recalculated, excluding 
the discarded experiments, with the standard error of 
the mean equal to 0.01mm. The resulting values D vary 
a little with a change of h, because of statistical fluctua- 
tions. 

To compare the experimental data, we therefore used 
dimensionless variables with /i-dependent units - the 



FIG. 4: Experimental results: The average drop destabilization 
length versus the distance from the wall, both normalized by 
the initial diameter of the drop. 

It is clear that the destabilization time and length of a 
drop are shorter for smaller distance h from the wall. In 
particular, the differences between values attained at h = 
2 mm and h = 10 mm are significantly larger than both 
statistical and systematic errors. In Sec. IVIi the results 
of the measurements will be extensively discussed and 
compared with results from computations based on the 
point-particle model. 
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FIG. 5: Experimental results: The average drop destabilization 
time T/r versus the distance from the wall h/D. Statistical er- 
rors (black) and systematic errors (gray) are indicated. 



IV. MODELING THE INTERNAL STRUCTURE OF 
SUSPENSION DROPS 



In the experiments, the suspension of particles is poly- 
disperse, as shown in Fig. [2j and the average number Nq 
of particles in a drop is greater than 7000. Numerical 
computation of the particle dynamics during a very long 
time, taking into account the size distribution, and for a 
large number of initial random configurations of parti- 
cles, located at different distances from the wall, would 
be very time consuming and complicated. Our goal is 
to concentrate on essential features of the drop evolu- 
tion, therefore we simplify the system. At the beginning 
of the experiments, a suspension drop is spherical with 
diameter D, it consists of iVo particles of diameters di, 
with i = 1, No, and an average of 25/im, and the drop 
is subject to the total gravitational force (weight minus 
buoyancy) F tot . 

In the numerical simulations, we use a simplified 
model of the drop internal structure. The average ini- 
tial diameters of a model drop and the real one are the 
same, D = 1.13mm, and the total gravitational forces 
are equal, F to t- However, a model drop consists of iden- 
tical particles of diameter d m . Moreover, to simplify the 
computations, the initial number of particles inside the 
model drop, N, is one order of magnitude smaller than 
inside the real one, Nq. 

To compensate for this difference, there are different 
possible choices of the particle diameter d m . In model 
I, we require that the volume fraction <fi m of the model 
drop is the same as the volume fraction <po of the real 
drops, 



or equivalently, 



[ ^ I 



(8) 



In model II, we assume equal mean diameters of par- 
ticles inside the model and the real drop, 



Model II: d„ 



, No 

— Yd, 

l—l 



(9) 



Both models are schematically shown in Figs. [6] and [7] 
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FIG. 6: Model I of the drop internal structure: equal volume 
fractions. 
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FIG. 7: Model II of the drop internal structure: equal mean 
particle diameters. 

Other values of d m , between those following from the 
model I and II, are also reasonable. The point is that the 
dynamics of the model drop only weakly depends on the 
choice of the particle diameter d m : as in model I, II or 
in between. To illustrate this property, we calculate the 
Stokes velocity of a single particle in a model drop, 



Vh 



3wr]d ri 



(10) 



Model I: 



(7) Then, we compare it to the initial drop velocity v, given 



5 



by Eq. ©. With the use of Eqs. © and ©, 



Vo = 5 
u 6Nd m /D 



0.02 for model I, 
0.05 for model II. 



(11) 



Therefore, in both models, and also for intermediate val- 
ues of d m , the particle Stokes velocity Vo is very small in 
comparison to the initial settling velocity of the drop v. 
A specific choice of d m is not relevant. 



V. THEORETICAL RESULTS 

In this paper, the theoretical approach is based on the 
point-particle model. In this way the evaluation of the 
particle dynamics is much simpler numerically and has 
a smaller number of parameters than, e.g., in case of the 
accurate multipole algorithm for spherical particles. 19 
The main advantage is that the point-particle dynamics 
does not depend on the particle radius, if it is carried out 
in the frame of reference moving with the Stokes veloc- 
ity^. 



A. Point-particle model close to a plane solid wall 

In a system with Reynolds number much smaller than 
unity (as in our experiments), the fluid velocity u(r) and 
pressure p(r) can be determined from the Stokes equa- 
tions. The fluid flow generated by motion of particles, 
which are subject to external non-hydrodynamic forces 
(e.g. gravitational ones), is often described within the 
point-particle modeliii^ In this case, the Stokes equa- 
tions take the form 



?/V 2 u(r) — Vp(r) 



V • u(r) = 



(12) 
(13) 



where r Q and F a are the position of a particle a, and the 
external force it exerts on the fluid, respectively. Here N 
is the number of particles, and r} is the dynamic viscosity 
of the fluid. The solution is the sum of the Green ten- 
sors appropriate for the specific geometry of the system 
and the boundary conditions. For an unbounded fluid, 
these conditions state that |u(r)| — > when |r| — > oo. In 
addition, for a fluid bounded by a solid wall located at 
z = 0, the stick-boundary condition has to be satisfied at 
the wall, i.e. |u(r)| = if r = [x, y, 0). 

We assume that all the point forces are identical, 
F Q = F, and parallel to the solid wall. We choose a sys- 
tem of coordinates in which F = (— F, 0, 0), with F > 0, 
as illustrated in Fig. [8] 

We choose the frame of reference moving with the 
Stokes velocity of a single point-particle in an un- 
bounded fluid. In this way, the dynamics of particles 
and the fluid flow are independent of the particle radius. 



image P' 




FIG. 8: The system, the coordinates and the notation. 

When the wall is present, the point-particle motion is de- 
termined by the following expressions for the particle ve- 
locities, a = 1, N, 

N 

v a =^T(f a ,r^F|T(r m |F. (14) 

Here the Blake tensor 20 

T(r a ,tp) = To(r Q(3 ) +T(r Qj8 /) (15) 

is the sum of the Oseen tensor (the Green tensor in the 
unbounded fluid), 



To(r a ^) 



1 



8irr]r al 3 



1 + 



and the tensor T(r Q , ( g/ ), defined as 

T(r a p>) • F = -T (r a p') F 
-2/^F-P- V r 



r a/3 _ 



3 ,T (r Q 



2/i 



+ 



87T?7 



F P V r 




(16) 



(17) 
(18) 

(19) 



The Oseen tensor, To(r Q( g), describes the interaction of a 
particle a with the particle j3 in an unbounded fluid, and 
the tensor T(r a p' ) describes the interaction of a particle a 
with the mirror image f3' of the particle /3. Here P = 1 
2zz is the reflection operator, hp is the distance between 
the particle /? and the wall, I the unit tensor, z the unit 
vector perpendicular to the wall, r Q7 = r a — r 7 , and r = 
|r|. The tensor T(r) accounts for the difference between 
the Green tensor for the fluid bounded by the wall and 
unbounded. 
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In Eq. O/ which defines the velocity of a point- 
particle a close to a wall, the first term is the sum of the 
fluid velocity fields generated by all the other particles at 
the position r a where the particle a is located. This term 
describes advection of a particle by the fluid flow gen- 
erated by all the other particles. The second term is the 
self contribution: it specifies the velocity of a single par- 
ticle interacting with the wall, in the frame of reference 
moving with the Stokes velocity of a single point-particle 
in an unbounded fluid. Since the wall is parallel to the 
force F, it follows that 



T(w) • F 



3F 



32h a TTT] 



(20) 



1,...,N is 



The dynamics of the point-particles a 
governed by the following system of first order ODEs 



dig 

dt 



(21) 



with the dependence of v Q on the positions rp of all par- 
ticles /3 = 1, N, given by Eq. (fl4)> . 



B. Evolution of a suspension drop 

In this section, we apply the point-particle model to 
describe the evolution of a suspension drop, with the 
same fluid outside and inside. N identical point forces 
F are distributed at random with a uniform TV-particle 
probability distribution inside a spherical volume of di- 
ameter D. In this paper, N = 700. Thirty such initial 
random configurations have been generated. 

The particle dynamics was determined based on the 
point-particle model close to a vertical wall, described in 
the previous section. The model was implemented nu- 
merically in our MATLAB code, using the variable order 
Adams-Bashforth-Moulton solver odell3 with the rela- 
tive accuracy of 0.1% and the absolute error of 10~ 6 . In 
the numerical simulations, the same set of initial random 
configurations was evolved for ten different distances 
h/D of the drop center from the wall (see Table [TTT1 for 
the list of values), and for a drop in an unbounded fluid. 

We introduce dimensionless quantities by taking D as 
the length unit, and v, given in Eq. |@), with 



Ftot = NF, 



(22) 



as the velocity unit, with F = |F|. The corresponding 
time unit r is defined by Eqs. 10 and ((22). It is known 
from the literature 2,21 that, in an unbounded fluid of 
viscosity r\, the settling velocity of a suspension drop 
of diameter D, subject to N identical, randomly dis- 
tributed point-forces F, is well approximated by the ex- 
pressions (ID and ll22l , if N is sufficiently large. For arbi- 
trary N, the exact theoretical expressions were obtained 
in Ref . |H|] by statistical averaging. An isolated suspen- 
sion drop settles with the same velocity v as a fluid drop 
of the same size and the same excess weight. 



In the simulation, we evaluated the dynamics of the 
particles in each suspension drop until Tk = 3920r. The 
generic pattern of the drop evolution was the same as 
in an unbounded fluid 3,5 and in our experiments. The 
statistics is shown in Table |lll] Initially, the suspension 
drops were spherical, with the particles randomly dis- 
tributed inside their volume. Later, the drops gradually 
lost single particles, which were left behind in a thin tail. 
During the motion, the drops expanded horizontally 
and contracted vertically, in a form of a horizontal torus. 
All the time, the particles recirculated inside the drop, 
moving faster in the inner parts of the drop, slower in 
the outer parts, coming back to the inner parts again, and 
so on. Suddenly, the torus bended, and broke into two 
(or sometimes three) fragments, which became spheri- 
cal and repeated the generic evolution pattern. Such a 
scenario happened for 77% of the total 330 suspension 
drops simulated numerically. In 98% of the destabiliza- 
tion events, the drop broke into two fragments; only 5 
drops out of 254 broke into three pieces. 23% of the 
drops did not break up before Tk = 3920r. Some of the 
particles were gradually lost from these drops, one by 
one, while the other particles inside the drops recircu- 
lated all the time. In Ref. JD], it was shown that in un- 
bounded fluid, the percentage of particle clouds which 
break up increases with the increase of the initial num- 
ber N of point-particles inside the cloud. 



TABLE III: The number of suspension drops with the indicated 
final stage of the evolution at T k = 3920r (the point-particle 
model). 



h/D 


destabilization 
into 2 droplets 


destabilization 
into 3 droplets 


no break-up 


total 


00 


24 





6 


30 


70 


16 


1 


13 


30 


30 


22 


1 


7 


30 


10 


24 





6 


30 


6.5 


24 





6 


30 


5 


22 





8 


30 


3.5 


24 


1 


5 


30 


2.5 


26 


1 


3 


30 


1.5 


19 


1 


10 


30 


1 


22 





8 


30 


0.75 


26 





4 


30 



In the following, we will discuss only the evolution 
of the drops which broke up. Analyzing the results of 
the numerical simulations, we use the same definition 
of the destabilization time T as in the experiments, see 
Sec. IIII Bl Within the point-particle model, the destabi- 
lization length L — VqT corresponds to the distance the 
drop traveled during time T, in the frame of reference 
moving with the Stokes velocity Vq of a single particle in 
an unbounded fluid. 

A typical statistics of destabilization lengths and times 
is shown in Fig. |9j where the initial distance of the drop 
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center from the wall h/D = 5. From Fig. [9^ and b it 
is clear that destabilization lengths, (L — VqT)/D, and 
times, T/t, of individual drops differ significantly from 
one suspension drop to the other. T/t can even change 
as much as almost one order of magnitude. Fluctuations 
of (L — VqT)/D are smaller, but still the values can differ 
even by a factor of 3. This effect is related to the chaotic 
nature of many-particle dynamics^ 
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FIG. 9: Destabilization lengths (a) and times (b,c) of individual 
drops, for h/D = 5. 



ages are taken over values corresponding to suspension 
drops which destabilized during the numerical compu- 
tation, i.e. at T < TJ. = 3920r, separately for each dis- 
tance h from the wall. Both destabilization length and 
time are smaller if the drop is closer to the wall. 
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FIG. 10: The average drop destabilization length (L — VqT)/D 
versus the distance from the wall h/D. The horizontal solid 
line and dashed lines correspond to the average destabilization 
length and the standard error of the mean, for an unbounded 
fluid. 



Another striking feature is the non-symmetric distri- 
bution of destabilization times, which more frequently 
attain smaller values, see Fig. 03. This asymmetry is re- 
duced if we look at the statistics of the natural logarithm 
ln(T/r) rather than T/t, as shown in Fig.|9}:. Therefore, 
the average destabilization time is evaluated as the ex- 
ponential function of the mean ln(T/r) throughout this 
paper (both in the numerical simulations and the exper- 
iments). 

The average numerical values of the destabilization 
length and time are listed in Table IIV1 and plotted in 
Figs. [TOWTTl Again, for simplicity, we use the same sym- 
bols L and T there to denote the mean values. The aver- 

TABLE IV: The average destabilization length, (L - V T)/D, 
and time, T/t, following from the point-particle model, for dif- 
ferent distances from the drop center to the wall, h/D. 



h/D 


(L - V Q T)/D 


T/t 
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185 ± 27 


601 ± 53 


70 


184 ± 24 


580 ± 63 


30 


189 ± 24 


600 ± 62 
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207 ± 24 
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160 ±21 
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FIG. 11: The average drop destabilization time T/t versus the 
distance from the wall h/D. The horizontal solid line and 
dashed lines correspond to the average value of T/t and stan- 
dard error of the mean, for an unbounded fluid. 



VI. COMPARISON OF THE EXPERIMENTS WITH THE 
POINT-PARTICLE MODEL 

A typical evolution of 8 suspension drops is shown in 
Fig. HZl For each drop, time is normalized by the desta- 
bilization time of this specific drop. In these units, the 
time separations between the subsequent frames 1-4 are 
equal to 0.28, and between the subsequent frames 4-6 are 
smaller, equal to 0.06, to show more precisely the drop 
just before and after the break-up. All the frames have 
the same size 4D x 4D. The gravitational field is vertical. 
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FIG. 12: Shape evolution of suspension drops settling in a fluid close to a vertical wall (left, A-D), and very far from a wall (right, 
E-H). A and B: Experiments at h/D — 1.75, with destabilization times T/t = 220 and T/t = 634, respectively. C and D: Point- 
particle simulations at h/D = 1.5, with T/r = 321 and T/t = 966, respectively. E and F: Experiments at h/D ss 88, with 
destabilization times T/t = 478 and T/r = 667, respectively. G and H: Point-particle simulations in unbounded fluid, with 
T/t = 274 and T/t = 1211, respectively. 



There are 4 pairs of drops shown in subsequent 
columns, A-B, C-D, E-F, and G-H. In each pair, the drops 
differ only by the initial configuration of the particles. 
The two drops in the pairs were selected such that the 
first one breaks up at a significantly shorter time T/t 
than the second one. The evolution patterns of the drops 
from the same pair are practically the same, if observed 
in time t normalized by the destabilization time T. 

In Fig. [121 we compare 4 experiments (highlighted in 
gray) with the corresponding 4 numerical simulations 
based on the point-particle model (white). The shape 
evolution observed experimentally is well approximated 
by the corresponding shape evaluated from the point- 
particle model. The shape evolution is quite similar for 
drops close to a vertical wall and far from it. In the ex- 
periments, the drop just before break-up is a bit more flat 
and wide when the distance from the wall is larger. 

The statistics of destabilization lengths and times for 



individual drops is shown in Fig. [13] In the experiment 
and in the model, there exist certain minimal values of 
both L and T below which the break-up does not occur, 
because drops need some time to change their shape ac- 
cording to the pattern shown in Fig.[12l in a similar way 
as it has been observed in point-particle simulations for 
an unbounded fluids 

The destabilization times tend to concentrate around 
values slightly above this limit, but some of them are sig- 
nificantly larger, even an order of magnitude. For the 
point particles, the range of values attained by both T 
and L is wider at larger distances h from the wall. No 
such tendency is observed in the experiments, where in 
general the spreading of values is smaller than in the 
simulations, slightly for L and significantly for T. The 
smaller spreading can be explained by the stabilizing ef- 
fect of the lubrication interaction between close surfaces 
of the particles. 23 "Evaporation" of particles from the 



9 



Point-particle model 




(L-V Q T)/D 



400 



1.2 



0.8 



Q 



0.6- 



0.4- 



0.2 



Experiment 



100 



200 
L/D 



300 400 



Point-particle model 



Experiment 



1.2 



0.8- 



0.6 



0.4 



0.2 



500 1000 1500 2000 2500 
T/x 




500 1000 1500 2000 2500 
T/t 



FIG. 13: Destabilization lengths and times of individual drops. Left: the point-particle model; right: experiments. 
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drop by statistical fluctuations is slowed down owing to 
this effect. In both the experiments and the simulations, 
the range of the relative differences between destabiliza- 
tion times of individual drops is very large, and it in- 
creases when h is getting smaller. The spreading seems 
to represent chaotic nature of the particle dynamics^ 

In Fig. [Mj we investigate the relation between the 
destabilization time and length, plotting T/t (in a loga- 
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FIG. 14: A universal linear scaling of \n(T/r) with L/D for 
all the destabilization events; top: the experiment, bottom: the 
point-particle model. [Color on-line. Top: h/D=1.75 (red), 4.63 
(cyan), 8.62 (green), 17.9 (blue), 26.8 (black), 43.1 (red-brown), 
87.7 (cherry). Bottom: h/D=0.75 (red), 1.0 (cyan), 1.5 (green), 
2.5 (blue), 3.5 (black), 5.0 (red-brown), 6.5 (cherry), 10 (ma- 
genta), 30 (battleship gray), 70 (dark pink), oo (purple).] 

rithmic scale) versus the destabilization length (in a lin- 
ear scale) for all individual suspension drops, separately 
for the experiment (top panel) and the point particle 
model (bottom panel). In each single panel, ln(T/r) as 
a function of the destabilization length is well approxi- 
mated by one linear relation (a least-squares fit) for all 
the drops. Deviations are observed only for the clos- 
est distances from the wall (h/D = 1.75 for the exper- 
iment, red on-line, and h/D = 0.75,1 for the point- 



particles, red and cyan on-line). The solid lines in Fig.lMl 
correspond to the least squares fits In T/t = (0.0030 ± 
0.0001)T/T> + (5.76 ± 0.03) for the experiments, and 
In T/t = (0.0079 ± 0.0002)(T - V T)/D + (5.06 ± 0.03) 
for the point-particle model. 

The experimental slope is less steep than the numer- 
ical one, and the constant term is higher. To compare 
the plots, the point-particle simulations need to be trans- 
formed from the frame of reference moving with the 
Stokes velocity Vq to the laboratory frame. This issue will 
be discussed at the end of this section. 

The experimental and numerical average destabiliza- 
tion lengths are compared in Fig. [15] In both cases, L/D 
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FIG. 15: The average destabilization length versus the inverse 
distance from the wall. Top: experiments, bottom: point- 
particle model. 

is approximately a linear function of the inverse distance 
from the wall, D/h. In detail, the least squares fits in 
Fig.HSareL/D = (-247±8)D/M-(304±2) (experiment, 
solid line), and (L - V T) /D = (-78 ± 9)D/h + (190 ± 5) 
(point-particles, dash-dotted line). For a given D/h, the 
experimental values of L/D decay faster with the de- 
creasing h than the numerical ones. The fast decay of 
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L/D can be related to the increased rate of the particle 
loss at smaller h, observed in the experiments. 

The average destabilization times, T/r, are compared 
in Fig. [16] The linear fits from Figs. [14] and [15] result in 
the curves plotted in Fig. [16] The experimental times are 
slightly larger than the numerical ones, but the differ- 
ence is smaller than the error bars. We can conclude that 
the experimental and numerical destabilization times 
are approximately the same, so that the influence of the 
drop distance from the wall on the destabilization time 
found in the experiments is captured well by the point- 
particle model. 
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FIG. 16: The average destabilization time versus the inverse 
distance from the wall. Gray: point-particle model, black: ex- 
periments. Dash-dotted and solid lines follow from fits shown 
in Figs.H4land[T5l 

The values of the destabilization times T/r, listed in 
Table [IV] and shown in Fig. [16] allow to estimate the 
shift, V T/D = (V /2v)(T/t), which transforms the 
point-particle positions at the moment of break-up, from 
the moving to the laboratory frame of reference. For 
the model of the drop internal structure, described in 
Sec. IIVI the parameter Vq/2v varies between 0.01 and 
0.03, depending on a specific choice of the particle di- 
ameter d m . The shift VqT/D to the laboratory frame of 
reference results in the destabilization length value L in- 
creased by around 3% and 8%, respectively This is a 
small change, and in the laboratory frame of reference, 
there is still a positive difference between the experimen- 
tal and numerical destabilization lengths. This effect will 
be discussed in the next section. 



VII. CONCLUSIONS 

In this paper, we have analyzed experimentally and 
numerically the influence of a vertical hard wall on the 
dynamics of particles settling under gravity in a viscous 
fluid at low-Reynolds-number. Initially, a large num- 
ber of close particles was randomly distributed inside 



a spherical volume. We have shown that the wall sig- 
nificantly affects the dynamics of this group of particles 
(called a suspension drop). 

The essential effect is that, in the presence of the wall, 
a drop breaks up faster and it travels a smaller distance 
as a cohesive entity. It has been shown that the average 
destabilization time T and distance L are smaller for a 
closer distance h of the drop center from the wall, with 
approximately linear dependence of T and L on D/h, for 
h larger or comparable to the drop diameter D. Other- 
wise, the typical pattern of the drop evolution is essen- 
tially the same with and without a wall, both in the ex- 
periments and the numerical simulations, as illustrated 
in Fig. El 

The investigations of this work were challenging, both 
experimentally and numerically. The measurements in 
glycerol are difficult, owing to the temperature depen- 
dence of the dynamic viscosity, and the tendency of glyc- 
erol to capture water from the air. In future investiga- 
tions, another fluid, e.g. a silicon oil, might be more prac- 
tical. 

The numerical results have been evaluated with the 
use of the point-particle (point-force) model, in the 
Stokes flow regime. Such a model has been widely ap- 
plied for an unbounded fluid. We have constructed and 
applied the point-particle model close to a vertical wall. 
For identical point forces, this model is very simple, be- 
cause the particle dynamics does not depend on the par- 
ticle radius. 

A challenge is also the comparison between the point- 
particle model and the experiment. It is known from the 
theoretical foundations and examples studied in the lit- 
erature 2 * that the model is too simple to accurately repro- 
duce the drop settling velocity. Therefore, the present 
paper can be also considered as a case study to estimate 
the accuracy of the point-particle model. 

For a given distance from the wall, the average desta- 
bilization time observed in the experiments and com- 
puted in the point-particle simulations is approximately 
the same. In the experiments, the destabilization length 
is greater than in the simulations, with the difference not 
exceeding 30% far from the walls, and systematically de- 
creasing to zero with decreasing distance between the 
drop and the wall. 

It seems that, in the experiments, the time-dependent 
drop velocities are higher than in the point-particle sim- 
ulations, although initially they were equal to each other 
- faster drops reach a longer distance, if they break at 
the same time. Drop velocities decrease with time, be- 
cause drops gradually loose particles. It seems that, in 
the experiments, a smaller fraction of particles is left be- 
hind the drop than in the simulations. In other words, it 
seems that the experimental drops tend to keep a larger 
fraction of particles, measured with respect to the initial 
number, and this can be the reason why they move faster 
than the numerical drops. This effect can be caused by 
hydrodynamic interactions between close particle sur- 
faces, which bind them together more strongly than in 
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the case of point-particles. Such a stabilizing effect of lu- 
brication between the particle surfaces has been shown 
and discussed e.g. in Ref. |23|l . 

To check this hypothesis, measurements of the num- 
ber of particles lost in the tail would be worthwhile, in 
comparison to the results of the point-particle simula- 
tions. The evolution of drops, which contain spherical 
particles, will be the subject of a future numerical study, 
based on the multipole algorithm of solving the Stokes 
equations and the numerical code hydromultipole^ 

An interesting result is also a statistical distribution of 
destabilization times and distances. Both T and L were 
shown to differ significantly from each other, even an or- 
der of magnitude, for different random configurations of 
the particles, owing to the chaotic nature of the particle 
dynamics. The smaller spreading of values in the exper- 
iments can be also explained by a stabilizing effect of the 
particle surfaces. 

The results obtained in this paper can be used in 
practical applications, such as sedimentation of micro- 
particles clouds, small milk or ink drops, in micro- 
channels or close to container walls. 



of L and T given in Table [TIJ was calculated from the 
standard formula, 



S x = 



n(n — 1) ^— ' 

. ' i=i 



(Al) 



During the experiments we observed large variations 
of the measured quantities from day to day. They might 
be related to small temperature variations, temperature 
gradients, convection of the fluid, water at the interface 
of the glycerol, or other reasons. Such a systematic error 
was estimated in the following way. First, we computed 
separately average values measured on a specific day k = 
1, ...,m, 



Xk 



1 lk 

T~ / Xi,k j 

2 — 1 



(A2) 



where i = 1, ik labels experiments performed on the 
day k. The corresponding standard deviation, 
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Appendix A: Estimation of statistical and systematic errors 
in the experiments 
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(A3) 



was then used to determine the weight Wk = 1 l^xk °^ 
the measurements performed on the day k, and to deter- 
mine, for a given value of h, the average value of \ and 
the standard deviation, 
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(A4) 
(A5) 



For a sequence of experiments i = l,...,n with fixed 
values of 0, D, and h, the standard error of the mean, 
X = - Ysi=i Xit (with x standing for L or T, and values 



Finally, we estimated the statistical error as the average 
of cr S y over all such distances h, for which measurements 
were taken on different days. 
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